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ABSTRACT 



Recently, the PAMELA, FERMI, HESS, and ATIC instruments have discov- 
ered interesting spectral features in the positron to total electron ratio and in 
the total electronic component of cosmic rays at high energy. These observations 
are studied with analytical models for the propagation of cosmic rays emanating 
from a discrete set of sources or from a smooth distribution of sources extending 
over the Galaxy. The contrast between the theoretically expected spectra from 
these two source distributions is seen to play a crucial role in the interpretation 
of the recent findings. It is shown that the positron to electron ratio, observed by 
PAMELA, may be fit over the entire energy range in the nested leaky-box model 
for cosmic ray propagation and that the ratio is expected to reach an asymptotic 
value of ~ 0.6 at high energies in all conventional models of cosmic ray propaga- 
tion. We also derive the spectral shape of the electrons and positrons expected 
from the annihilation of dark matter in the Galaxy and show that the spectral 
shape of the peak will provide important information, not only regarding the 
mass of the dark matter particles, but also regarding their spatial distribution. 
The spectrum of secondary positrons and electrons, calculated with the nested 
leaky-box model, is subtracted from the spectrum of the electronic component 
to determine the spectrum of primary electrons emerging from the cosmic ray 
sources. This spectrum is analyzed in terms of contributions from a set of dis- 
crete sources sprinkled across the Galactic volume. Our analysis suggests the 
possible presence of electrons, accelerated in high Mach number shocks, in the 
cosmic rays. 
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1. Introduction 

Direct observation of the cosmic ray electronic component (which includes both 
electrons and positrons) without any charge discrimination dates back to the early 1960s, 
and since that time, the energy range and the sensitivity of the observations have increased 
systematically. To date, we have at hand da ta from three new instruments, FERMI 



( lAbdo et al. 2009h . ATIC fiChang et al. 20081 ) . and HESS (Aharonian et al. 2008, 2009) 



that have the requisite sensitivity to measure, with good statistical accuracy, the spectrum 
of the total electronic component {e~ + e~^) well into the TeV region. The reported spectrum 
in the region ii^ > 10 GeV may be parameterized as 

FtiE) = UE-^^ (1) 

with Ft = 3.05 up to ~ 1 TeV, and the value of Ft increases to ~ 3.9 at higher energies. 

Furthermore, the ATIC collabortion has reported a narrow enhancement of the 
intensities around 600 GeV, superimposed on the smooth spectrum. We reproduce in fig. 
1 their observations and a compilation of the results of other measurements. The spectral 
slope below 10 GeV progressively flattens to a slope of ~1.7. Uncertainties in the estimates 
of the interstellar fluxes are introduced due to solar modulation effects below an energy of 
a few GeV. We also show, in the same figure, a smooth fit to all the data that we adopt for 
some of the calculations. 

The electronic component in cosmic rays, because of its interactions with radiation 
fields such as starlight, the microwave background, and magnetic fields in the Galaxy, has 



been particular 



y useful in underst anding the origins and propagation of energetic particles 



in the Galaxy (ICowsik et al. 19661 ). This and other early considerations of the effects on 
the spectral shape of the cosmic ray electrons were carried out in the context of a smooth 
distribution of cosmic ray sources in the Galaxy, and the transport was described within 
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the framework of the leaky-box model (ICowsik et al. 19671 ). In its initial formulation, the 



cosmic ray residence time in the Galaxy was assumed to be independent of energy for ^ 1 
GeV. However, improved measurements showed that the ratio of secondary nuclei to their 
parent primary nuclei decreased with increasing energy beyond ~ 1 GeV. This observation 
naturally gave rise to the suggestion that the residence time of cosmic rays in the Galaxy 
could indeed be energy dependent and decrease with increasing energy. According to this 
model, the anisotropy of the cosmic ray fluxes will correspondingly increase with energy 
and will exceed the observational bounds at very high energies. Motivated by the need 
to alleviate this difficulty and by general considerations, the nested leaky-box model was 
proposed, which took into consideration the effects of storage of cosmic rays in a small 
bubble or cocoon surrounding the compact sources of cosmic rays (Cowsik & Wilson 1973, 
1975). Eventhough the residence time of cosmic rays in the cocoon was energy dependent 
to generate the energy dependent part of the secondary to primary ratio of cosmic ray 
nuclei, subsequent to their leakage from the cocoon, the residence time of cosmic rays in the 
Galaxy was assumed to be independant of energy. Thus, the anisotropy remained sensibly 
constant up to very high energies, ~ 10^ — 10® GeV. We provide some details of this model 
later in this paper. 

The discrete nature of the sources of cosmic rays was addressed within the context 
of a diffusion mode 



and the implication to the electronic component was worked out 
(ICowsik fc Lee 19791 ) to show that we need sources within a few hundred parsecs in order 
to reproduce the observed spectrum of electrons up to ~ 1 TeV. Here, the sources of cosmic 
rays were distributed randomly in a thin disk, and the diffusion region was a thick disk. 
The boundary condition was that the density of cosmic rays vanished at the upper and 
lower surfaces of the thick disk, signifying free escape. The need for the proximity of cosmic 
ray sources to explain the electron spectrum at high energies was reiterated by similar 
analyses by Nishimura and others (Nishimura et al. 1997; Atoyan et al. 1995). In many of 
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these papers, the complexities introduced by the boundary conditions on the two planes 
bounding the diffusion volume is avoided by assuming that the diffusion volume is of infinite 
thickness. 

In this paper, we would also like to keep the Green's function simpler than that 
provided by Cowsik and Lee (1979), but we would like to still have the option of leakage 
of the particles from the Galaxy. To this end, we include a term signifying simple 
leakage similar to that in the leaky-box model. The considerable success of this model in 
reproducing the energy spectra from the cosmic ray sources and the fact that the Green's 
function integrates to the leaky-box model for a uniform distribution of sources surrounding 
the observation point justifies this approach. Apart from the transparency with which the 
results reported are derived here, we also note that many of the results are indeed generic 
to any reasonable model of cosmic ray propagation. 

It may be noted that this model fits the cosmic ray observations such as the ratio of 
secondary to primary nuclei in the cosmic rays, the positron fraction, the upper bounds on 
anisotropy, and the total electron spectrum. Whereas we refer our readers to the papers 
cited here (Cowsik et al. 1967; Cowsik & Wilson 1973, 1975) for the full details of the 
model, we provide here all the details needed for the discussion of the recent data on the 
electronic component. We show that the spectral slope of the electronic component at the 
highest energies will be dominated by secondaries generated by the nuclear component of 
cosmic rays through collisions with the interstellar matter. This is because as long as the 
cosmic ray sources are discrete, the electrons directly accelerated by them would be cut off 
sharply due to synchrotron radiation and the inverse Compton effect at energies determined 
by the distance to the nearest source. 



The recent observations of positrons in cosmic rays by PAMELA (lAdriani et al. 2009al ) 



has created much excitement because of the possible interpretation of these observations as 
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annihilation of dark matter in the Galaxy. Their observations and the results from earlier 
measurements are compiled in fig. 2. Even at the outset, it is worth pointing out that the 
positron fraction observed by PAMELA has a higher statistical significance and a similar 
energy dependence in comparison with the earlier measurements, but falls markedly lower 
over the entire energy band of the measurements. It appears difficult to ascribe the lower 
ratios observed by PAMELA to the effects of solar modulation, as they are seen even up 
to ~ 20 GeV. The positron fraction at ~ 1.64 GeV was measured by PAMELA to be 
~ 0.0673, which decreases to ~ 0.0483 at ~ 6.8 GeV and thereupon increases monotonically, 
reaching a value of 0.137 at a mean energy of 82.55 GeV. It is this monotonic increase 
that is being called anomalous, as it does not conform to the prediction of the leaky-box 
model of the cosmic ray pro pagation with energy dependent residence time in the Galaxy 



( iMoskalenko fc Strong 19981 ). Accordingly, this paper begins with a review of the cosmic 
ray propagation models. The theoretical expectation that the positron fraction should level 
off or gently decrease with energy is specific to the assumption that the ratio of secondary 
to primary nuclei, like B/C, will continue to fall up to very high energies. In interpreting 
the positron fraction, it is important to note that there is a fundamental difference between 
kinematics of the generation of positrons and secondary nuclei by the interaction of primary 
cosmic ray nuclei; whereas the secondary nuclei emerge from the collisions with nearly the 
same energy per nucleon as the primary nuclei that generate them, the positron produced 
through the interactions of cosmic ray nuclei carries away only a few percent of the energy 
per nucleon of the primary. 

This fact leads to the very interesting consequence that the secondary positron 
spectrum in the interstellar space has the same spectral slope as the primary nuclei in 
the context of the nested leaky-box model even though the ratio of secondary to primary 
nuclei, such as B/C, follows the observed decrease with energy. Of central relevence to the 
discussion of the positron fraction is the spectrum of the cosmic ray nucleonic component 
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that is dominated by protons with some neutrons coming in bound as He and other nuclei. 
The nucleon spectrum may be represented as 



Fr,{E) = I^E- 



(2) 



where r„ ^ 2.6 — 2.7 (We adopt r„ = 2.65.) in the energy region 1 GeV/nucleon to 
~ 10® GeV/nucleon, beyond which the slope may increase to ~3. The p/n ratio effectively 
determines the e^/e~ ratio generated by cosmic rays via interactions with matter in the 
sources and the interstellar medium through which they propagate before they leak out of 



the Galaxy. The theoretical calcu 



ation of /e 



generated through nuclear interactions of 



cosmic ray nuclei yield ^ 1.5 — 2 (jProtheroe 19821 ). On the other hand, the observations of 



the ratio (IHayakawa 19691 ) in cosmic rays gives fi~^/fi ~ 1.3. 



In recent efforts to explain the positron fraction and the "bump" in the total electron 
spectrum observed by ATIC around 600 GeV, it has been s uggested that astrophysical 



objects such as a nearby gamma ray burst source (jloka 20081 ). one or more pulsars (Hooper 
et al. 2008a; Kawanaka et al. 2009; Malyshev et al. 2009; Pohl 2008; Profumo 2008; 
Yiiksel et al. 2008), or a few nearby supernova remnants (Dado & Dar 2009; Fujita et al. 
2009; Shaviv et al. 2009) may be the cause of the "anomaly". Most of the recent effort 
has been in looking for a dark matter explanation for these observations. Supersymmetric 
dark matter such as the neutralino may explain the excess in the observed total electron 
spectrum, but there may be a need for excessive clumping or an exceedingly local source 
(Hooper et al. 2008b; Ishi wata et al. 2009). Th is clumping can be avoided if one considers 
Kaluza-Klein dark matter jnooper et al. 2009l) . An interesting dark matter explanation 



comes from the introduction of new force carriers (lArkani-Hamed et al. 20091 ) which allows 
for a Somerfeld enhancement and boosts the annihilation rate to leptons. The new particles 
may even b e able to explain the annual modulation of the signal seen it the DAMA 



experiment (IBernabei et al. 20081 ). Various other dark matter explanations have been 
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formulated to explain th e recent findings, and we refer the reader to (jProfumo 20081 ) and 



( iMertsch &: Sarkar 20091 ) for more detailed lists of the recent studies. 



In section 2, we provide a brief overview of the models for cosmic ray propagation 
and in section 3, we discuss the PAMELA results, first in a model independent way and 
then compare it with the expectations of the various models. In section 4, we subtract 
the secondary positron and electron spectra from the spectrum of the total electronic 
component observed by FERMI, HESS, ATIC, and other experiments to obtain the 
spectrum of electrons arising exclusively due to input from various discrete sources of cosmic 
rays sprinkled over the Galaxy. We then analyze this spectrum to derive the distance to 
the nearest source and the mean spacing between the cosmic ray sources in the Galaxy. In 
section 5, we analyze the excess in the total electron spectrum by considering a 5-function 
input spectrum and the spectrum expected from shock acceleration. Finally, section 6 is 
devoted to a discussion of the main results of this paper and related matters. 

This p aper is an extension o f work presented at the 31st International Cosmic Ray 



Conferece fICowsik fc Burch 2009h . 



2. Brief Overview of Models of Cosmic Ray Propagation 

It is generally accepted that cosmic rays are accelerated in a large number of discrete 
sources distributed in the Galaxy, and the cosmic rays propagate from these sources moving 
along randomly oriented trajectories, akin to diffusion, until they leak away from the 
Galaxy. During such a propagation, the cosmic rays might interact with matter surrounding 
the sources, interstellar matter, radiation, and magnetic fields. Any secondaries generated 
through such interactions, if charged, will be confined by the interstellar magnetic fields 
and will therefore follow the same kind of random paths as the primaries before escaping 
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from the Galaxy. The various propagation models are characterized by the specific form 



chosen for the "vacuum path length distribution" (ICowsik et al. 19671 ) which describes 
the probability P{t) that the cosmic rays spend in any given region, such as a cocoon 
surrounding the sources, or in the general interstellar medium before escaping into the 
intergalactic space. The term "vacuum" emphasizes the fact that in specifying P{t), one 
considers hypothetical particles which do not suffer interactions or lose energy during 
propagation. The effects of these processes are to be computed based on the distribution of 
paths P{t). 



2.1. The Leaky-Box Model 



In its simplest original form (jCowsik et al. 19671 ). one assumes that P{t) has a broad 



distribution with significant amplitude near t = 0, exemplified by a simple exponential 
function 

P{t) = e"*/" (3) 

where r is called the escape lifetime of the cosmic rays. In the original version, r was 
assumed to be sensibly independent of energy beyond ~ 1 — 2 GeV. Thereafter, since the 
discovery that the ratio of the fiuxes of secondary cosmic rays to those of the primaries was 
a decreasing function of energy, r was considered to decrease with energy to accommodate 
the observations. We summarize in fig. 3 and fig. 4 the available observations. The crucial 
issue here is how r behaves at energies beyond 10 — 20 GeV where the observations have 
low statistical significance or are non-existent at much higher energies. Most conventional 
models today (Moskalenko & Strong 1998; Strong et al. 2007) assume that 

Model A : ta{E) ~ toE"^, E>2 GeV/n (4) 

with A ^ 0.4 — 0.5. Such an extrapolation to high energies is shown as a dotted line in 
fig. 3 and fig. 4. To understand the spectra of the nuclear component in this model let 
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the rate at which sources inject cosmic ray nuclei, say carbon, into the galactic volume, per 
unit volume, unit energy interval, and unit time be given by 

S,{E) = K,E-^^ (5) 

Assuming that the leakage from the Galaxy dominates their loss, the equilibrium 
spectrum in the interstellar space is given by 

F,{E) = S,{E)ta{E) = K,T,E-^^^+''^ = I,E-^^ (6) 

In order to match the observed spectrum of carbon nuclei with ~ 2.65, we should choose 
/3c = Tc — A 2.15, i.e. a relatively flat injection spectrum. 

Now, secondary nuclei, like B, are generated at the same energy per nucleon as their 
progenitor C, at a rate 

Sb{E) = cNh(Jb,c{E)F,{E) (7) 

where Nh is the density of target nucleon in the interstellar medium. The B production 
cross section aB,c{E) is nearly independent of energy and 

Sb{E) = KbE-^^ (8) 

where Kb — cNh(Tb,c{E)KcTq — cNh(Tb,cIc- The equilibrium sepctrum of B is given by 

Fb{E) = Sb{E)t^E-^ = XBro£;-(^=+^) = cNHaB,cKy^E-^^^+^'^'^ (9) 

Comparing this expression for Fb{E) with that for Fc{E) in equation 6, one sees that 
Fb/Fc ratio falls as E'^. In other words, the ratio of secondary nuclei to their primaries 
just follows the energy dependence of the leakage lifetime in the leaky-box model. It is 
worth noting that in this model the anisotropics in the cosmic ray fluxes increase as E^, 
becoming unacceptably large for E ~ 10^ GeV. 
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In an attempt to overcome this difficulty we may assume that the residence time 
follows the energy dependence of the secondary to primary ratio, say up to ~ 10 — 20 GeV 
and levels off theeafter to a constant value. That is, we may assume that t{E) becomes 
nearly constant at high energies 



where tq is a constant. This is shown as a solid line in fig. 3 and fig. 4. The constancy 
of T ~ tg at high energies predicts a lower level of anisotropy of cosmic rays without any 
increase with energy. However, in this model the spectral slope of the primary nuclei, 
r„ = /3„ + A, for < 10 GeV, and r„ = for > 10 GeV. Thus in order to match the 
constant value of r„ ^ 2.65, we will have to assume that /3„ = 2.65 — A a; 2.15 for E < 10 
GeV, and /3„ = 2.65 for E :$> 10 GeV. 

We will refer to the two models described briefly here as leaky-box model A and 
leaky-box model B respectively. 



An alternate way of accommodating the falling secondary to primary ratio is in the 
context of the nested leaky-box model (Cowsik & Wilson 1973, 1975). Here, one assumes 
that subsequent to acceleration, cosmic rays at lower energies spend some time in a 
cocoon-like region surrounding the sources, interacting with the matter there to generate 
some of the secondaries. The residence time in the source region is assumed to be energy 
dependent, decreasing with increasing energy. On the other hand, once these cosmic rays 
enter the general interstellar medium, their subsequent transport becomes independent of 




Eb - 10 GeV 



(10) 



2.2. The Nested Leaky-Box Model 
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energy and the residence time becomes equal to tq- 

{Ts(E) ~ tb(E) - To 1 GeV < E < 10 GeV 

sy J By J G ^^^^ 
Ts{E) ~ To 10 GeV <E<10^ GeV 

The net effect of the interactions in these two regions is to generate the correct ratio of the 
fluxes of secondary nuclei to those of their primaries, as shown in fig. 3 and fig. 4 marked 
with labels C. 

Let the rate at which a cosmic ray source accelerates primary nuclei and injects them 
into a cocoon surrounding the source be represented as Sn{E) and given by 

SniE) = q^E-^- (12) 

The total number of cosmic rays within the cocoon is determined by the balance between 
the injection rate and the leakage characterized by the hfetime Ts{E). 

UE) = Sn{E)Ts{E) = q^E---Ts{E). (13) 

The nuclear component leaks out of the cocoon at a rate that is inversely proportional to 
the residence time Ts{E). Thus each source injects cosmic ray nuclei into the interstellar 
medium at a rate given by 

fn{E)/T,{E) = Sn{E) = QnE-'^- (14) 

If the spatial number density of the cosmic ray sources is represented by u, then the 
injection rate per unit volume is given by 

Sr,{E) = uqnE-'^-. (15) 

Note that this spectral exponent q,„ is identical to that of the input into the cocoons. In 
the above analysis, we have neglected the slight decrease in the value of qn due to nuclear 
interactions with the mateiral in the cocoon; more exact expressions are available elsewhere 
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(Cowsik and Wilson 1973, 1975), which are also apphcable to the rest of the discussions in 
this section. Now, the leakage hfetime from the Galaxy tg, in the nested leaky-box model, 
is independent of energy, and the equilibrium spectrum of primary nuclei Fn{E) in the 
Galaxy is given by 

Fn{E)^Sn{E)TG^iyqnrGE-"- = InE-'^-, i.e. = r„ (16) 

Note that this represents the spectral form intrinsic to the acceleration process. In constrast 
with this result for primary nuclei, note that the rate of generation of secondary nuclei 
within the cocoon, such as B, by the spallation of primary nuclei, such as C, is given by 

SBiE) = cnHaB,cfc{E) = cnH<JB,cqcrsiE)E-''^ (17) 

This leads to a spectrum of secondary nuclei inside the cocoon fsiE) given by 

fB(E) = SB(E)rs{E) = cnHaB,cqcrl(E)E--^ (18) 

and the effective rate of injection of secondary nuclei by the cocoons into the insterstellar 
medium becomes 

SBiiE) = vUE)/ts{E) = cnHaB,ci^qcTs{E)E-'''^. (19) 

This rate of injection is just due to the spallation taking place in the sources and reflects the 
energy dependence oiTs{E) superimposed on the primary spectrum ~ E~°'''. The spallation 
of primary nuclei in collision with interstellar matter also generates secondary nuclei at the 
rate 

Sb2{E) = cNHaB,cFc{E) = c7V^(7s,ci^gcrG^""= (20) 

These two injection processes for secondaries will lead to an equilibrium spectrum of 
secondary nuclei, Fb{E), 

Fb{E) = {Sbi{E) + Sb2{E)}tg = caB,cqci^TG{nHTs{E) + NhTg}E-^''. (21) 
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As noted before, these remarks are made with the assumption that the production 
of the tertiary nuclei and higher order effects are small. However, it is straightforward to 
derive the exact expression in a more general case (Cowsik & Wilson 1973, 1975). Note 
that the expression {nHT{E) + NhTs} in equation 21 has been chosen so as to reproduce 
the observed B/C ratio in cosmic rays. 

It is important to note that in model C, the spectrum of the primary nuclei as observed 
by various experiments has the same spectral index as that generated by the acceleration 
process, both being ~ —2.65. The absence of significant anisotropy of cosmic rays at high 
energies is also easily understood. The residence time of cosmic rays in the Galaxy is 
essentially constant up to ~ 100 TeV. This will predict that anisotropy will not increase 
significantly with energy at least up to this energy. 

2.3. Spectrum of Secondary Electrons and Positrons in Cosmic Ray 

Propagation Models 

The generation of electrons and positrons in the interactions of the cosmic ray nuclear 
component occurs through the production of mesons, mainly pions, which decay to muons, 
which in turn decay into electrons or positrons, transferring, on the average, a fraction of 
about 0.05 of the energy per nucleon of the primary. This is in contrast with the production 
of secondary nuclei, such as boron from the collision of carbon nuclei, where boron emerges 
with almost the same energy per nucleon as the primary carbon nucleus. This difference in 
their production characteristics leads to nearly identical source spectra Sn- and iS'„+ for the 
secondary electrons and positrons ~ E~^" in all the three models (A, B, and C). To see 
this, we note that the rate of generation of secondary positrons and electrons at any energy 
E is proportional to the flux density of the nucleonic component at £"„ ~ E/0.05 ~ 20E and 
the collision rate. Thus in all the models (A, B, and C) the rate of production of positrons 
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and electrons in the interstellar medium is given by 

Sn± cx F„(20E) ~ E-^- (22) 

In considering the nested leaky-box model C, we should in principle add the generation 
of secondaries inside the cocoons. This rate is proportional to 

Sn±c{E) oc WOE) ~ g„E-""r5(20^) ^0 E :^ 1 GeV (23) 

In writing equation 23, we have used equation 13 to define /„ and note from fig. 3 and fig. 
4 that Tg becomes very small compared with tq for En = 20E > 20 GeV. Thus the injection 
spectra of positions and electrons for E > 1 GeV are essentially the same in all the three 
models. 

On the other hand, their equilibrium spectra Fn+{E) and Fn-{E) are markedly 
different in the three models. For energies below ~ 100 GeV, where the energy losses due to 
synchrotron radiation and inverse Compton scattering on radiation fields are not important, 
the three models generate the spectra noted below: 

Model A: ~ 5„+(E)rA(^) ~ roS"(^"+^^ (24) 

Model B : F„+ ~ Sn+{E)TB{E) ~ ro^"(^"+^\ E < 10 GeV 

~ S„+{E)tg ~ E-^", E>10 GeV 

Model G: Fn+ ~ S'„+(E)rG ~ E > 1 GeV 

The spectra for the electrons are similar to those given in equation 24, except that 
because of the dominance of the protons in the cosmic ray beam, the production rate of 
electrons is somewhat lower, with 



This ratio r; is theoretically estimated 



Tom the characteristics of high energy interactions 



to be ~0.5 flMoskalenko fc Strong 19981 ): on the other hand, the direct observation of the 
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fi / ratio indicates a value of ~ 0.8 (IHayakawa 19691 ) . In either case, r] is essentially 



independent of energy beyond a few GeV. 

Thus we see that, at high energies {E ^ 10 GeV), in the leaky-box model B and in the 
nested leaky-box model C, the secondary positron and electron spectra are power laws with 
indices Tf, = F^, i.e. equal to that of the spectrum of the nuclear component in cosmic rays. 
At very high energies, the energy losses due to synchrotron emission and inverse Compton 
scattering will steepen these spectra to F„+ = E~^^'=~^^\ We will discuss this aspect further 
in section 4. 



3. Analysis of the Positron Fraction Observed by PAMELA 

It is useful to write the observed positron fraction R{E) in terms of the various 
components: 

Here, F„+ and F^- represent the positron and electron spectra generated as secondaries of 
the nuclear component of cosmic rays, and F^- is the spectrum of electrons resulting from 
direct acceleration in the sources. Note that there is no direct acceleration of positrons in 
the source. It is convenient sometimes to work with the inverse of R{E) given by 

This allows one to find the spectrum of electrons generated by the sources Fg- as 

F,-{E) = [P{E) - (1 + v)]K+{E). (28) 

We show in fig. 5 the net secondary spectrum Fn±{E) = (1 + n)Fn+{E) calculated in model 
C, along with the spectrum of electrons resulting from direct acceleration in the sources, 
Fe-{E), which is obtained by simply substracting Fn±{E) from each of the observed data 
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points displayed in fig. 1. We also show in the same figure a solid line representing F^-^E) 
calculated using equation 28 and a smooth fit to the PAMELA observations; this line and 
the data points in fig. 5 represent the spectrum of electrons generated exclusively through 
acceleration in various cosmic ray sources in the Galaxy. 

Alternatively, we may just assume the functional form for Fn+{E) given by the various 
propagation models and calculate the positron fraction by dividing Fn+{E) by Ft{E), the 
total spectrum of electrons measured by FERMI, HESS, and other experiments, 

R.(E) ^ (29) 

which is shown in fig. [6] along with the data from PAMELA and other experiments. The 
normalization of the theoretical curves is such as to provide the best possible fit to the 
three models A, B, and C described earlier. (This normalization may indeed be explicitly 
calculated as it is proportional to ta{E), tb{E), and tq respectively for the three models 
and depends on the density of matter in the propagation region, the spectral fiux of the 
nuclear component, the cross section for meson production, decay kinematics, etc.) In 
depicting the three model curves, we have included the effect of the radiative energy losses 
at high energies. 

Comparison of the theoretical expectations of the positron fraction with the PAMELA 
data indicates that model A provides a r ather poor fit to the observations, as already noted 



by several authors flAdriani et al. 2009a] ). A careful calculation of the positron fraction 
under the general assumptions of model A was carried out over a decade ago by Moskalenko 
and Strong (1998). Our estimates here are essentially the same as that derived by them. 
Even though both model B and model C predict nearly identical injection spectra, the 
equilibrium spectra at low energies in these two models differ drastically with each other. 
On the other hand, both these models predict identical equilibrium spectra at high energies 
for E ^ 10 GeV. The good fit to observed the positron fraction shows that the residence 
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time of cosmic rays is essentially independent of energy for E > 10 GeV, a constancy that 
is expected to continue up to ~ 10^ GeV. 

In choosing between model B and model C, the latter is preferred from considerations 
of the spectra of primary nuclei as well. This is because for a simple power law input from 
the sources having a form Sn{E) ~ E~^, model B would be expected to yield a spectral 
form Fn{E) ~ E'^^^^^ below ~ 10 GeV and E'^ at higher energies. On the other hand, 
the observed spectra of all the primary nucleonic components are simple power laws of 
slope ~ E^^'^^ with no changes of slope in the tens of GeV region. Thus we conclude that 
the nested leaky-box model provides good fit with the PAMELA data and is preferred also 
from consideration of other cosmic ray observations. 

In order to derive the behavior of the positron fraction at very high energies, beyond 
the PAMELA domain, we need to study the possible form of the spectrum of primary 
electrons generated by the cosmic ray sources, Fe^i{E). The key point to keep in mind 
here is that there exists a discrete set of compact cosmic ray sources, which are distributed 
randomly over the Galaxy. This discrete nature of the sources and its impact on the form 
of the electron spectrum from a source F^^i^E) have been discussed earlier, notably by 
Cowsik and Lee (1979) and by Nishimura et al. (1997). The diffuse flux from a transient 
source at a distance r peaks at a time t ^ r^/Qn, where k is the diffusion constant. This 
finite t ime of propagation i nduces a sharp cutoff in Fe^i{E) at high energies due to radiative 



losses fl Cowsik fc Lee 19771 ). On the other hand, the source functions 5'„+ or S'„-, which 
correspond to the product of cosmic ray flux and the density of matter in the interstellar 
medium, have smooth spatial distributions, and accordingly, lead to a simple steepening 
of the spectrum to the form ~ E~'^l^^^\ Thus, at the highest energies, the secondary 
component Fn+{E) and Fn-{E) will dominate over Ff,i{E), the spectrum of electrons 
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generated by the sources, so that 

R(E) ~ ^ /f^? , , = ^— ~ 0.6. (30) 

Thus we expect, on very general considerations, R{E) to increase at high energies and reach 
a plateau at a level of (1 + n)~^ ~ 0.6, which is dictated by the p to n ratio in the interstellar 
cosmic ray flux and by the nature of nuclear interactions at very high energies. Based on 
the discussions in the section 4, we expect the positron ratio to reach the asymptotic value 
at > 5 TeV. 



4. The Spectrum of Primciry Electrons Generated by Cosmic Ray Sources 

For the purposes of this discussion, we write down the transport equation for cosmic 

rays as 

!?^_,V=iV + ^ = (31) 
at T 

where N is the number density of cosmic ray particles at a distance r from the source at 
time t, K is the diffusion constant, r = tg is the mean lifetime for the escape of particles 
from the Galaxy (taken to be independent of energy in the nested leaky-box model), and 
Q is the source term. Setting Q — 5{t)5{r), we get the Green's function G{r,t) for the 
transport 

G{T,t) = {A7:Kt)-'/^exp[ - ^ - (32) 

Note that in writing equation 31, we have not included the term representing the energy 
loss for the electrons V • {hE'^N). This is because the energy loss due to synchrotron 
and Compton processes takes away the energy of the electrons in small steps so that the 
energy loss may be treated as continuous without any stocasticity. This loss is customarily 
described by the equation 

^ = -hE^ (33) 
dt ^ ' 
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or 



^ E{t = 0) - E{t) ^ Eq-E 

bE{t^O)E{t) bEoE ^ ' 

where Eq is the energy E{t — 0) and E = E{t). Thus, with the subsidiary condition given in 
equation 33, the effect of energy loss may be completely taken into account (Cowsik & Lee 
1979). We note also that for a uniform distribution of sources surrounding the observation 
point, this Green's function leads to the leaky-box model, with an exponential "path length 
distribution" 

P{t) = J G{r,t)4Trr^dr ^ e~*/\ (35) 

Consider now a source located at a distance ri, which accelerates particles continuously 
to a spectrum of the form 

Qi{E) = Q,E-f' forE<E,. (36) 
The spectrum that will be observed is given by 

I bE bEx r, r, cs 

Fa{E, n) = / QeE-l'il - hEtf-^G{n, t)dt. (37) 

The conceptual meaning of equation 37 is clear: equations 33 and 34 imply that the energy 
of the electron or positron at the instant of production, E{t = 0), is given by 

^<' = "' = T^ 

and a unit bandwidth in energy at E{t — 0) gets compressed to a bandwidth (1 — bEty, 
enhancing the flux density per unit energy interval by this factor. The upper limit of 
t = {1/bE — 1/bEx) coincides with E{t = 0) = E^., the maximum energy up to which the 
sources accelerate cosmic ray electrons. For the present purposes, we may take E^ ^ oo. 

Before we proceed further, we should fix the value of b to be used. To this end, we note 
that for E ^ 100 GeV, the scattering cross section of the starlight by the electrons becomes 
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small, as described by the Kline-Nishina formula. Thus we need to consider the effects of 
the scattering only on the microwave background at 2.7 K and synchrotron radiation in the 
magnetic fields in the thick disk, which is the region of cosmic ray confinement. The value 
of b is given by the formula 

b = 3.22 X 10-^{Wph in eV/cm-^) + 7.9 x IO'^HIq^^^^ GeV'^Myr-^ 

^ 8.5 X 10-^ + 7.1 X 10-^ 1.56 x 10~^GeV-^Myr-\ (39) 

We have taken Wph ~ 0.25ey cm~^ and < H >= S/iG. 

We start with a discussion of the spectra for a smooth distribution of sources such as 
that for the secondary positrons and electrons in model C; their spectral shape is given by 

pl/bE 

Fs{E) = / K±E-^{l-bEtf-'^e-*'^dt (40) 
Jo 

r 1 

~ K±tE-^ for E<^Ec^ — 

br 

~ K±b-^E-^^+^^ for E:^E^. 

It is this form of that we adopted in the theoretical estimates for the secondary positron 
spectrum and positron fraction in section 2. 

The contribution of secondary electrons and positrons Fn± to the observed spectrum 
of the total electronic component Ft{E) in cosmic rays is simply 

F„±(E) = F^+{E) + Fr,-{E) = (1 + r?)F„+(E) (41) 

where Fn+{E) has the form Fs{E) given by equation 40 and is normahzed to yield the 
correct positron fraction as described in section 2. In order to obtain the spectrum of 
electrons generated exclusively by the cosmic ray sources, we need to subtract Fn±{E) from 
the total observed spectrum 

F,-{E)^Ft{E)-F„±{E). (42) 
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We show in fig. 5 Fg-{E) obtained by substracting F„±(£') from the data from ATIC, 
HESS, FERMI, and other experiments. 

In order to understand the nature and the distribution of the sources of cosmic ray 
electrons, we first show in fig. 7 the spectrum Fd{E, ri) for various values of ri and a power 
law input with F = 3 calculated using equation 37. We note that these spectra are similar 
to those calculated by Cowsik and Lee (1979), and they appear to fit the observations for 
ri ~ 100 pc. However, the spectrum Fg-(E) has been generated by the inputs from an 
ensemble of discrete cosmic ray sources in the Galaxy located at distances rj. Assuming 
that the sources are randomly distributed with a mean spacing equal to the distance to the 
nearest source, the average distance to the i*'* source is 

< Tj >fs nVi (43) 

where ri is the distance to the nearest source. The theoretical expectation for the 
equilibrium spectrum generated by the ensemble of cosmic ray sources for a power law 
input is given by 

N N i/hE , 

Fo{E) = V F,{E, n) = T ^{1 - hEtf^-^G{< n >, t)dt (44) 

i=i i=i •^*=o ^ ' 

where N is the total number of sources which contribute to the local cosmic ray density. 

The spectrum Fd{E) has been calculated for various values of ri, and the results are 
displayed in fig. 8 along with the spectrum Fg-{E) estimated in equation 42 from the 
observations of FERMI, HESS, PAMELA and other experiments. We note that for the 
distance to the nearest source of ~ 100 — 200 pc, and for a typical separation of ~ 100 — 200 
pc between the sources, the theoretical calculations match the observations reasonably well, 
except that they fall short of the intensities beyond a few hundred GeV, leaving behind a 
narrow spectral feature centered around 600 GeV. This result, derived here more carefully 
and with better data, confirms the earlier analysis of Cowsik and Lee (1979) and Nishimura 
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et al. (1997) that the spectrum of the cosmic ray electronic component at high energies 
shows that there are cosmic ray sources at distances of the order of ~ lOOpc from the solar 
system and that the typical spacing between them is expected to be the same. 



5. Narrow Spectral Features in the PrimEiry Electron Spectrum 

The narrow spectral feature may be isolated by subtracting from the spectrum of the 
primary electrons Fg-{E), the expected spectrum Fd{E) due to power law inputs and 
diffusion from these sources: 

n,{E) = F,-{E) - Fr){E) (45) 

This is displayed in fig. 9, and this excess or more specifically the sharp spectral features 
observed by ATIC, has been ascribed in fMo or in part to products of dark matter 
annihilation, as mentioned in the introduction. In this section, we discuss this excess in 
terms of two input spectra, one a (5-function in energy and the other a fiat spectrum ~ 
such as that expected for acceleration at planar shocks of high Mach number. 



5.1. 5-function Input Spectrum 

The injection spectrum is assumed to be of the form 

Q(t = 0,E„) = 5(E(t = 0)-^„), (46) 

which could in principle represent the spectrum arising from annihilation of dark matter. 
Without any leakage from the Galaxy, such a spectrum will evolve as 

Note that the above spectrum integrates over E to unity for all values of t. Now, for a 
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continuous injection and an energy independent leakage, the spectrum is given by 

F.,(£,£J = ^e.p-(^). (48) 

We display in fig. 10 the spectrum calculated for Ea — 1200 GeV for a spatially smooth 
distribution of sources and compare it with the narrow feature ne{E) estimated in equation 
45. In fig. 11, we provide further examples of the theoretical spectra for various values 
of Ea. Note that is positive for E < = ^ and negative for E > Here, 
Ec — ^ 600 GeV. Thus, for a (^-function input at Ea > 300 GeV, there is no peak in 
the observed spectrum at the energy value E — Ea- Therefore, the narrow feature seen at 
~ 600 GeV is unlikely to be generated by an input 5-function source spatially extended 
and smoothly distributed about the solar system. Incidentally note that if the value of b is 
about a factor of two smaller than our estimate of 1.56 x 10"^ GeV~^Myr~^, we can get a 
broad peak at ~ 600 GeV, but then the secondary component will extending into the HESS 
region as shown in fig. 5 (dotted line). 

Next we consider a 5-function input from single discrete source located at various 
distances rj. This leads to an observed spectrum 

E"^ / hE E \3/2 / hr'^EE E —E\ 

We display in fig 12, several examples of Fsd for different values of r^. This spectrum 
displays a peak at 

For small r^, this peak will be sharp near Ea, but with increasing r^, this peak shifts to lower 
energies and becomes broader. As before, we sum over the sources at various distances to 
find their net contribution Fsd{E) and show this in fig. 13 for different values of the mean 
spacing between the sources. Such a composite spectrum is quite broad, with a very slight 
maximum around 200 GeV and does not provide a good fit to the feature at 600 GeV. 
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5.2. Shock Acceleration 

Acceleration of particles at planar shocks of high Mach number generate spectra which 
are power laws of index ~ 2 

QshockiE) ~ QoE-^ for E < E,. (51) 

The upper cutoff energy E^ and the precise shape of the cutoff depends on the nature 
of the shock, radiative energy losses, and other such features. Such an input leads to an 
equilibrium spectrum 

F,{E)^^{l-e-^y (52) 

This is shown in fig. 14, and it appears to reproduce ne{E), the difference between the 
primary electron spectrum Fe_ {E) and that expected by a set of discrete cosmic ray sources 
FdiE). 

6. Discussion and Related Matters 

The main result that emerges from the present analysis is that the nested leaky-box 
model provides a good fit to the positron fraction observed by PAMELA. The model is also 
consistent with other observations of cosmic rays. Until good measurements of the positron 
fraction were available, there was no easy means of choosing amongst various models. The 
fact that the nuclear secondaries, such as Li, Be, and B, emerge from nuclear interactions 
with essentially the same energy per nucleon as their parents, C, N, and O, was the main 
cause for this degeneracy. However, the fact that the positrons carry, on the average, a 
fraction of only about 0.05 of the energy per nucleon of their nuclear primaries breaks 
this degeneracy, allowing a choice to be made. Improvements of the measurements of the 
spectra of both the secondary nuclei and of positrons will help in fixing, more firmly, the 
parameters of the nested leaky-box model. Note that this model predicts that the p/p ratio 



-26- 



( lAdriani et al. 2009bl ) should nearly be constant at ii^ ^ 10 GeV. Also, the gamma rays 
generated by cosmic rays inside the cocoons through the 7r° 27 process will have, on the 
average, a steep spectrum E~^Ts{E) ~ i^-C^n+A) qx > \ GeV. In closing, we note that 
as the centenary of the discovery of cosmic rays is approaching us, the field of cosmic ray 
studies is becoming more vibrant and is intimately getting connected with wider aspects of 
research in astrophysics. 

The recent measurements of the electronic component of cosmic rays by the PAMELA, 
HESS, and FERMI groups has opened up many interesting new avenues to discuss cosmic 
ray propagation. Some of these were discussed or hinted at in this paper. Further d iscussions 



will include an anysis of the p/p ratio as seen by PAMELA (lAdriani et al. 2009bl ). 



Acknolwedgement:It is a pleasure to acknowledge the detailed discussions with 
Professor Martin A. Israel and the hospitality at the Raman Research Institute where part 
of this manscript was drafted. 



-27- 



1000 



^. 100 - 



LL 
LU 



10 - 



Smooth Fit 
HESS Fit 
PAMEIJ^ Fit 



^ , — ^ 












* 


AMS -L 




m 


ATIC 







BETS 




* 


HEAT 




A 


HEAT94 




T 


HESS 




* 


Kobayashi 






PPB-BETS 




o 


FERMI 






HESS {LOW E) 


10 




100 

Energy [GeVJ 



1000 



Fig. 1. — The compilation of measurements of the total electron spectrum Ft{E), the red 
dashed line represents the total spectrum of the electronic component calculated using the 
positron ratio measured by PAMELA. The blue dashed line is a fit to the HESS data, and 
the solid line is a smooth fit to the total electron spectrum. 
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Fig. 2. — The positron fraction measured by PAMELA along with the earher measurements 
are shown. Note that PAMELA data has a similar energy dependence, but is markedly lower 
than the earlier measurements. 
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Fi g. 3. — The observed B/C secondary to primary ratio is plotted [points from a compilation 
in (jStrong et al. 20071 )] along with the power law extrapolation at high energies (dotted line, 
Model A), a constant extrapolation (solid line, Model B), and a two-component fit (dot- 
dashed hne followed by solid hne for E > 50 GeV, Model C). 
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Fig. 4. — The o 



jserved (Sc+Ti+V )/Fe secondary to primary ratio is plotted [points from a 



compilation in (IStrong et al. 20071 )] along with the power law extrapolation at high energies 
(dotted line, Model A), a constant extrapolation (solid line, Model B), and a two-component 
fit (dot-dashed fine followed by solid fine for E > 50 GeV, Model C). 
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Fig. 5. — We display here the primary spectrum of cosmic ray electrons Fe^{E) generated ex- 
clusively through acceleration of electrons in the cosmic ray sources obtained by substracting 
the secondary positrons and electron (solid line) from the measurements of the total elec- 
tronic component Ft{E) by HESS, ATIC, FERMI and other experiments. The sum of the 
secondary component and F^{E) the primary component will add up to Ft{E) showin in 
fig. 1 
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Fig. 6. — The theoretically calculated positron fraction for models A, B, and C are com- 
pared with observations. All calculations are normalized at ~ 10 GeV, to the PAMELA 
measurements. 
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Fig. 7. — The primary electron spectra due to a single source at various distances from the 
source with = 5 TeV compared to the primary electron spectrum, [ri = 0.1 kpc (solid 
line), ri = 0.2 kpc (dashed line), ri = 0.5 kpc (dotted line), ri = 1.0 kpc (dot-dashed line)]. 
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Fig. 8. — The theoretical primary electron spectra resulting from an ensemble cosmic ray 
sources for various values of the mean spacing and = 5 TeV is compared with the 
primary electron spectrum Ff.-{E). The mean spacing between the sources is taken to be 
< r >= 0.1 kpc (solid line), < r >= 0.2 kpc (dashed line), < r >= 0.5 kpc (dotted line), 
and < r >i= 1.0 kpc (dot-dashed line). 
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Fig. 9. — The excess primary electrons obtained after subtracting the expected contribution 
from discrete sources estimated with a mean spacing of ~ 0.1 kpc from the primary electron 
spectrum Ff.^{E). This is shown both as data points and as a smooth fit through the data. 
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Fig. 10. — The spectrum for a 5-function input from continuous distribution of sources with 
Ea = 1200 GeV . Note the peak at ~ 300 GeV^ Ec/2 is expected for all > 300 GeV (see 
fig. 11 
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Fig. 11. — Further examples of the spectra of electrons expected for 5-function inputs from 
a continuous spatial distribution of sources. 
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Fig. 12. — The spectrum of electrons expected for a (5-function input, diffusing spatially from 
a single source situated at various distances. 
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Fig. 13. — The spectrum of electrons arising due to a (5-function input from a discrete set of 
sources located at various distances with mean spacing as indicated, and calculated assuming 
diffusive transport. 
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Fig. 14. — The equilibrium spectrum of electrons arising from shocks to input spectra ~ E'"^ 
up to various cutoff energies between 5 — 10 TeV is compared with the narrow feature in the 
observed spectrum. 
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